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Verifying relationship between Height and Spacing, in Barchan Dunes 
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We have investigated the relationship between height and spacing of Barchan dunes 
which the coupled map lattice model numerically generates. There is a scaling relation 
between them and the values of the scaling exponents agree well with real dunes' values. 
The values of these scaling exponents are the same for both steady states and transient 
states. 
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1. Introduction 

There are many serious problems all over the world. One of them is how to control the behavior 
of dunes. So the problem of controlling dunes should be solved and it has recently begun to be 
researched quantitatively by a number of physicists. 

There is an experimental observation; height H and spacing L of dunes have a relation 

L ~ H a , (1) 

where a = 0.58 ~ 1.92. 1 ) Since it is a non-trivial observation, some theoretical analysis is necessary. 
Unfortunately there are no theoretical approaches for this relation. Instead we have a powerful 
method; numerical analysis using computer. In this research, we would like to investigate the 
validity of this relation by computer simulations, "Coupled Map Lattice (CML) model". It was 
developed by Kaneko, 2 ** and was applied to researches of dunes by Nishimori and Ouchi. 3 ' 4 ) 

A crescent-shaped barchan dune is frequently observed in a desert. We analyzed barchan dunes 
by the numerical simulation with CML model. In this paper, "dune" means "barchan dune". 

In the next section, we will discuss the relationship between height and spacing of real dunes. 
In the third section, we will show some results by computer simulation. Some discussions will be 
included in the fourth section. 
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2. Experimental Law and Assumption 



When there are N dunes, the nth dune is supposed to have height h n and mass m n . It is 
expected that they have similar shapes with each other, because granular matter has an angle of 
repose. However, if the sand particles are blown by the wind, smaller dune may lose more sand 
than the larger ones. They will break the similarity of shapes. Therefore, the mass of a dune is 
not proportional to 3rd power of height. Sauermann et al. 7 ^ show that 

m n = 7 h 2 n \ (2) 

where, 7 is a constant number, a total sand mass M tota i is given by 

N 

M total =jJ2 h n 4 - ( 3 ) 

n 

It is assumed that this Mtotai is shared by individual dunes. 

Barchan dunes are formed when M tota i is relatively small and the wind direction does not 
fluctuate. Thus we can choose two axes. One is parallel to the wind direction and the other 
is perpendicular to the same one. Along one of these axes, the distance between dunes can be 
measured as L11 or L±. Lm is a spacing of dunes along the direction parallel to the wind direction, 
and L± is a spacing of dunes along the direction perpendicular to the wind direction. Here, we 
know the area S n which nth dune occupies, 

S n ~L\\L L . (4) 

There is another expression of S n using the area Stotai occupied by N dunes. 

o _ Mtotai _ Stotai _ ft _ ^ 2 /, \2 4 /r\ 

where S to tai is an area of this system, I is size of the observed area. 
We assume that two relations, 

L\\ ~ (K) a , (6) 

L± ~ (h n f, (7) 

the occupied area by the nth dune is 

S n a L\\L ± oc (h n ) a (h r y = {h n ) a+ P. (8) 
(5) is compared with (8), the exponent of (8) is expected as 



a + (3 = 2.4. 



(9) 
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3. Numerical experiments and Results 

We used CML model which was arranged by Nishimori and Ouchi 3 ' 4 ) in order to apply to 
computer simulation. At first, we prepare a 1000 x 1000 lattice and give an uniform random 
number [0, X] to all sites. These random numbers correspond to the height for each site; 

h(i,j)e[0,x], (10) 

X is the parameter which controls the total amount of the sand, and h(i,j) is an amount of sand 
at a site (i, j). 

Saltation and creep are important when considering a dynamics of dune. Saltation is a flying 
of sand by the wind. Creep is that sand rolls and falls with gravity. Although they had better to 
be considered exactly, they are approximated for simplicity. 

Flight distance L s and the amount of the sand which flies by saltation q s have the relationship 
with Ah, which is a difference of height along the leeward direction, 6 ) 

q s = -tanh(A/i) + 1.3, (11) 

L s = tanh(Aft) + 1. (12) 

We consider that wind is blowing to the positive direction of the i axis, 

Ah(i,j) = h(i + l,j)-h(i,j). (13) 

Thus saltation is decided by only landform. 

Next, we consider creep, which is regarded as "diffusion of sand ". We introduce an arbitrary 
diffusion constant D, and the amount of diffusion of sand q c is described as 

q c = DAh(i±l,j±l), (14) 

using (13) and (14), we get 

/i(i±l,.7±l)-/i(i±l,.7±l) + g c , (15) 

h(i,j) -> h(i,j) - q c . (16) 

But, we must not forget "critical angle" of sand; it is an "angle of repose" of sand. According to 
the observation of sand in a hourglass, it is about 34°. The creep does not occur, if an inclination 
is equal to 34° or is less than that. Since tan 34° ~ 0.67, only when 

h(i±l,j±l)-h(i,j)>0.67 (17) 

creep occurs. 

We repeat the whole process. One step contains all of these. With this CML model, a de- 
scription of dynamics of dunes becomes much easier than a real situation. Nishimori and Ouchi 3 ' 4 ) 
have already reported that their model reproduced dune patterns qualitatively. We are interested 
in how similar it is quantitatively to the real dunes. 
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3.1 Scaling Relation in Steady States 

Dunes get to the steady state after some transient period. In order to realize this stationary 
state, we have iterated whole processes over sufficiently long period. In this case, it is 3,000 steps. 

We will discuss the average height of dunes. We do not employ simple averages over individual 
dunes. Instead we use the total height of sand, (H st ) which can be defined as 

(H st ) = H st /X = h st (i,j) 2 /X. (18) 

Since it is assumed that system size is fixed, (H st ) is proportional to the average of height (see 
Appendix). After some period, the system becomes steady states. h s t(i,j) expresses the quantity 
of the sand in a site at that time. Figure 1 shows the relation between X and a quantity 

(H s t) which is proportional to an average height of these steady dunes. There is a linear relation 




Fig. 1. Relation between mass of sand of average of height of saturated dunes 
between X and (H st ) as 

(H st ) /10 5 = (4.64 ± 0.07) A + (2.70 ± 0.07), (19) 

where CML model generates dune-like patterns only when this relation stands. Therefore this 
equation enables us to judge whether the system converges to the steady state or not. 

We examine both the quantity (H s t) which is proportional to height and the spacing^ L s t\\ 
which is along the direction parallel to the wind direction (see Fig.2), and get 

L st |, ~ (F st )°- 8±0 - 5 . (20) 

On the other hand, a quantity (H st ) has the relation with the spacing L st ± along the direction 
perpendicular to the wind direction, 

L st ± ~ (Hst) 1 ** - 1 . (21) 
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Fig. 2. Relation between an average height and a spacing of dunes (left;parallel to the wind direction, right;vertical 
to the wind direction) 

Substituting the exponents of (20) and (21) into the relation (8), we get 

S n oc <i^} 2 - 4±0 - 5 . (22) 
This result is consistent with the results obtained in real dunes as described in the section 2. 

3.2 Scaling Relation in Transient States 

Next we consider the scaling relation in transient states. Here, we employ two methods for this 
investigation. One is for a fixed time and another is time series scaling. 

3.2.1 One Fixed Time 

At first, we consider a certain fixed time in a transient state (t = 2, 000 steps), and investigate 
the relation between the spacing of dunes Lf x \\ which is along the direction parallel to wind direction 
and a quantity (Hf x ) (see Fig.3). A quantity (Hf x ) is defined as, 

(H fx )=Y,h fx {i,j) 2 /X. (23) 

Since system-size is fixed, (Hf x ) is proportional to the average of height at one fixed time in 
transient states. hf x (i,j) is the quantity of the sand in a site (i,j) at that time. And we get, 

L M ~ (H fx r±°\ (24) 

Next, we compute spacing Lf x ± which is along the direction perpendicular to the wind direc- 
tion, and get 

L fx ± ~ ///, (25) 
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Fig. 3. Relation an average height and spacing of dunes (left;parallel to the wind direction,right;vertical to the wind 
direction) 



Substituting the exponents of (24) and (25) into the relation (8), we get 

a + f3 = 2.3 ±0.3, 

therefore, 

S n oc //„. 

This result also agrees to the argument in Sec. 2. 



(26) 
(27) 



3.2.2 Time Series Scaling 

The transitional time period varies with the initial quantity of sand. During this period, state 
of dunes continues to change. Next, we analyze dunes in transient states, by time series scaling 
method. 

The total sand height in a certain time t in transient states is set to H tota i (t) . And the quantity 
H(t) which is proportional to a height of dunes in this states can be defined as 



Hit) = H total (t)/X = Y J h {i , j) {tf/X. 



(28) 



1,3 



Here, h^j^(t) expresses the height a site has in a certain time t. And A is a parameter which 
controls initial quantity of sand. We normalize variables H it) and L as 

t H(t) 



(t,H(t)) 



(H st ) si (H st ) 



), 



(29) 



-Lic+ J -'s 



(30) 

J st ^ st 

Where 8, 5' are scaling indices, and L st means a spacing of dunes in a steady states. We assume 
that a scaling relation between time t and height H(t) as follows, 



(H st ) J{ (H st )S>~ { (H st y 



(t « (H st ) S ), 



(31) 
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and also the same scaling relation between t and L(t) is assumed 

^ = /(4r) ~ (4)' (* « ^), (32) 

^s* ^st -^st 



removing t from (31) and (32), we get 



H{t) (H st ) L* 
L(t)7 L 7 { (H st ) 5 



^(TTT^Y- (33) 



J st 

The right hand side of (33) is constant. Therefore, the relation H(t) and L{t) is given by 

H(t)~L(t)? — ► L(t) ~ H{t) e T. (34) 
The parameter for the spacing along the parallel direction to the wind direction is 

a = (35) 
and for the spacing along the perpendicular direction to the wind direction is 

(3= € -f. (36) 




Fig. 4. Time series scaling of the height of dunes 

e ||( e _!_) is a parameter e' when a parallel (perpendicular) direction to the wind direction is con- 
sidered. In our simulations, three exponents (e, £||,e±) are (see Fig. 4 and Fig.5) 

e = 0.11 ±0.02 (37) 
ey = 0.08 ± 0.01 (38) 
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Fig. 5. Time series scaling of the spacing of dunes(left: parallel to the wind direction, right: vertical to the wind 
direction). L st || L st ± is a spacing along direction parallel to the wind direction in steady states. 



e_L = 0.16 ±0.01. (39) 

Here, we note that a ((3) is the exponent for parallel (perpendicular) direction to the wind direction 
and it is estimated as follows, 

eii 0.08 ±0.01 

a = 7 = im±ao2^ ' 7±0 ' 2 ' (40) 

tJi = 016*001 
' e 0.11 ±0.02 v ; 

Here we again give the value 

S n <x H(tf 2±0A . (42) 
Again, this value does not disagree with the values obtained previously. 



4. Discussions 

We have estimated the exponent a + (3 = 2.4 using several methods. This is the first estimation 
using numerical simulation. These values are consistent with the value by the argument developed 
from of Sauermann's observation 7 ) in Sec. 2, as shown in Table. 1. 

In the present paper, we have discussed a quantity which is proportional to an average height 
and spacing of dunes. All our results agree with Lancaster's observations. ^ The CML model 
approximates dynamics of dunes. Pay attention to the fact that we did not treat exactly relationship 
between sand and wind. But as you saw, our results agreed well with the real dunes. Thus, exact 
discussions about the relationship between wind and sand, for example a sand flux and a distance 
of flight, may not be so important. A sand flux and a distance of sand flight is designed by how 
wind blows over a dune, and how wind blows over a dune is designed by the landform. To tell the 
truth, even if we know about only landform, we can guess sand flux q s and distance of sand flight 
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state of dunes 


parameter 


value from observation 


2.4 


steady states 


2.4 ±0.5 


transient states (fixed time) 2.3 ± 0.3 


transient states (scaling) 2.2 ± 0.4 



Table I. Summary of this study 



L s . Similarly, creep was treated as a diffusion of sand. We did not consider repellent force and so 
on which must be considered when we deal with real sand creep. 

Thus, we can consider that exact theory about wind velocity, sand flux and so on is not so 
necessary. There seems to exist an universality class about the relationship among landform, sand 
flux and distance of sand flight. Thus the exact creep theory can be replaced with the diffusion 
theory. Now we suggest that CML model is suitable model when studying an average height and 
a spacing of dunes. 
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Appendix: Definition of average height 

Consider two differently-shaped dunes with the equal amount of sand. The following figures 
show the cross section along leeward direction. For simplicity, we assume a simple sinusoidal shape 
for each dune. The total height H sum and H' sum are 




H sum = I h(i)di = h n I sinaidi = — -, (AT) 



[ h(i)di 
Jo 

/ h(i)di = h' n / ° sina'idi = ^f. (A-2) 
Jo Jo a 

These integrated values are equal because the amount of sand is equal. And following relations are 
realized, 



and 

H' = 



2h n 2ti, 
7> ' 



(A-3) 



and 

h n : a = h' n : a'. (A-4) 

The average height of dunes (H) and (H') becomes equal to the integrated value divided by system 
size, 

(H) . ^, (A.5) 

and 

00 = §• (A-6) 

When (A-4) is considered, (H) = (H') is realized and it turns out that this calculation cannot tell 
us the height of dunes, h n or h' n . When the value of each site is squared, we get the integrations 

H s t = 



= / h(i) 2 di = h 2 n [ a sm 2 aidi = (A-7) 
Jo Jo 2a 



and 



H' st = J' = h' n 2 jj sin 2 a'idi = (A 



Since the average squared-height of dunes H st and H' st become equal to the integrated value divided 
by system size, next relations are realized, 

TQ=% (A.9) 
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and 2 

^=lf < A ' 10 > 

When (A-4) is considered, 

H7f.Hi t = h n :h' n (A-ll) 
is concluded. On the other hand, it is also that the ratio of these integration values is, 

H st : H' st = h n : h n >. (A- 12) 

From (A-7),(A-8), (A-ll) and (A-12), 

H7f.Hl t = H st :H' st = f h(i) 2 di : /' h(i)' 2 di. (A-13) 

Jo Jo 

Thus the ratio of average heights of dunes can be estimated by squared quantity of each site. 
If we define quantity (H st ) 

(Hst) = (A-14) 
and consider that total sand mass is the product of X and the system size £, 

H sum = XL (A-15) 

we get 

(H st ) = - h n t (A- 16) 
If i is fixed, we can use (H st ) as the quantity proportional to h n . 



